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The decay of a quasiparticle in an isolated quantum dot is considered. At relatively small time the 
probability to find the system in the initial state decays exponentially: P{t) ~ exp(— Ft), in accor- 
dance with the golden rule. However, the contributions to P(t) accounting for the discreteness of final 
three-particle states, five-particle states, etc. decay much slower being ~ (As/r)" exp(— rt/(2n-|- 1)) 
for 2n + 1 final particles. Here A3 <C F is the level spacing for three-particle states available via 
the direct decay. These corrections are dominant at large enough time and slow down the decay 
to become ln(P) ~ asymptotically. P{t) fluctuates strongly in this regime and the analytical 
formula for the distribution W{P) is found. 
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The problem of the decay of a quasiparticle in closed 
mesoscopic systems has attracted a permanent interest 
last years |l|-pi]|] . The very existence of the (irreversible) 
decay becomes nontrivial in the case of a discrete spec- 
trum. In a many particle system one should in addition 
understand the role of the whole hierarchy of the dis- 
creteness, ranging from the single-particle levels (level 
spacing A) to the exponentially dense spectrum of com- 
plicated many-particle excitations. It was first realized 
that the decay into two particles and one hole, relevant 
for an infinite geometry, does not apply for the low ex- 
citation energies in the finite system. Namely, for a dif- 
fusive quantum dot, the coupling to the three-particle 
states describes the broadening of the quasiparticle only 
for excitation energy e ^ \fg^, where 5 ^ 1 is the di- 
mensionless conductance Utilization of the ideas 

of delocalization in the Fock space for the quasiparticle 
decay in Rcf. 0] has stimulated further active research 
in this field 1^-^. However, the main interest was con- 
centrated around the limited range of excitation energies 
yjgl\n{g) < e/A < y/g. In spite of the beauty of the 
applied theoretical approaches experimentally it may be 
difficult to look for the effects sensitive to the logarithm 
of large parameter (see however Therefore the goal 

of this paper is to consider the effect of the discreteness 
of the spectrum for much higher energy, where the usual 
three-particle decay is expected to determine the quasi- 
particle lifetime. As we will see, even for e ^ y^A the 
decay in the confined geometry at large enough time is 
strongly modified and slowed down. 

The simple fermionic Hamiltonian models the effect of 
electron interaction in closed quantum dot (in the nu- 
clear physics an analogous Hamiltonian was introduced 
already a long time ago 0,n2|) 
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The single-particle energies Si are randomly distributed 
around e;^ = with mean level spacing A. The depth of 



the Fermi sea is always much larger than the excitation 
energy. The Gaussian random two-particle interaction 
has zero mean and variance given by IslT 
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The conductance 5 3> 1 measured in units of e^//i is a 
large parameter in our problem. We are interested in the 
stability of simple single-particle excitations above the 
ground state of H. Due to the weak interaction the low 
lying excited states are almost unperturbed. However, 
the level spacing for three-particle states (two particles 
and one hole) accessible for the direct decay decreases 
while increasing the excitation energy e [pi 
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and at e > y^A the matrix element of the interaction 
(^ exceeds A3. Now the Fermi golden rule is applied to 
find the width of the quasiparticle 
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The natural quantity characterising the time evolution 
of the single-particle state is the return probability P{t), 
which is the probability to find a system in the initial 
state |0) after time t 
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Here we have introduced Gqo (e) , the diagonal matrix ele- 
ment of the single-particle Greens function of the Hamil- 
tonian (|l|). 

The authors of Ref . [Q have used the Cayley tree model 
in order to describe the quasi-particle disintegration into 
many-particle states. One site of the lattice (tree) is as- 
sociated with the single-particle state which is connected 
with many three-particle states (sites) . Each of the three- 
particle states is connected with a number of five-particle 
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states, and so on. On the Cayley tree the exact Greens 
function is given by the simple formula M] 
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where e\ ^ e ^ i\ with some small A. Each new denom- 
inator of the continued fraction corresponds to a mixing 
with the next generation (a set of states in the Fock space 
having one more excited particle and one more hole). 



and so on. is the matrix element connecting the single- 
particle and three particle states and Vnm is the matrix 
element connecting the n-th three-particle and m-th five- 
particle states. All matrix elements satisfy Eq. (|^). Not 
all effects of the two-particle interaction are taken into 
account by the Cayley tree model Nevertheless, one 
may associate any term of the perturbative expansion 
(in V) of Eq. (||) with a certain class of true diagrams of 
many-body perturbation theory. The diagrams notpre- 
sent in Eq. (|^) should be taken into account separately, 
but they will not change the results of this paper. In 
particular, the interaction of the initial particle with the 
three-particle states is completely taken into account by 
the Cayley tree. Replacing the first sum in the denomi- 
nator in Eq. (0) by the integral over averages, one gets 



Goo(£) -eo(e) = (£-eo + *r/2)- 



where the width 



r = 2 
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essentially does not depend on the secondary width r„. 
The pole of Qq leads to the exponential decay of the re- 
turn probability ~ exp(— Ft) and it is hard to change 
this result while considering separately G{e) and G*{e'). 
The slow tail of -P(i) at large t may come only from terms 
in the correlation function Goo(e)GQQ(e') that are non- 
analytic in e — e'. 

Let at first stage forget about the decay of secondary 
states in Eq. (i.e. we put Vnm = 0). Now the ex- 
pansion of the Greens function (^ around the averaged 
value (Eq. (0)) gives 



Gqq — Gq + Gq 



dx 
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Again ex — £ + iA. Only the correlation of the second 
term {...} in Goo(£) with its analog in G'j5o(e') may lead 
to the correction singular in e — e'. It is convenient to 
replace the first sum in Eq. (^ also by an integral via 



J2 F{En) JJ2S{x- En)F{x)dx. Now the only non- 
trivial part in the calculation of the correlation function 
is the averaging 
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We suppose that £„ are independent random variables 
(Poisson statistics). In fact the problem of the decay of 
a simple state into the quasi-continuum of discrete states 
(just like our decay into three particles) has already been 
considered by several authors [p^^. In particular, the 
case of constant interaction and Wigner-Dyson statistics 
for secondary states is considered in . 

With the use of Eqs. (|,|l^) one finds {Gq = ^o(e')) 

Goo(e)GSo(£') = ^0^^* + (^0^^*)'^ , (11) 

e — e — 2iA A3 

which after substitution into Eq. (||) gives (for Gaussian 
random interaction V'^ — 3A^/2(7^) 
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Here 2A works effectively as the constant intrinsic width 
of the three-particle excited states. After the initial state 
spreads over the '--^ F/A3 directly connected states, the 
usual Breit-Wigner decay saturates. It is important for 
this result that the interaction matrix elements and(or) 
the three-particle interlevel spacings essentially fluctuate. 
Otherwise the coherent decay may proceed up to much 
larger times (A;^^) and even may become reversible. 

The many particle origin of the final states is in fact 
ignored in Eq. (|2|). In the real problem the energy and 
width of the "intermediate" final state are given by the 
sum of energies and widths of the constituent particles 
En = £1 + £2 + £3 « £ and F„ = - E £? §)■ In 

order to take into account this energy dependence of the 
secondary width one should return to Eq. (^ and con- 
sider the three-particle density of states as the sum over 
individual particle energies. This means that in Eq. (|lT 
the following replacement should be made 
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The return probability now takes the form 



p ^ g-rt ^ 27V3 A3 1 ^-rt/3 
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Here the states from the second generation are aUowed 
to decay and their widths are distributed within r/3 < 
r„ < r. However, the fluctuations of the widths, and 
even more important, correlations of the fluctuations, are 
stifl not taken into account. The second term in the r.h.s. 
of Eq. ( p^ ) is not smaU necessarily and for Ft > (3/2)L 
even becomes dominant. Here we introduced a special 
notation for the large logarithm L = ln(r/A3). How- 
ever even at these times the prefactor l/Tt in the return 
probability is only logarithmically small. We will neglect 
the contributions of this type in the calculation of further 
corrections to the return probability. 

In the way similar to the derivation of Eq. we may 
extract from the Greens function the term responsible 
for fluctuations in the coupling of three-particle to five- 
particle states 

Goo = go + - + g'oT.SnV^ (15) 

n 

Here A„ ~ A3 is the interval between the five-particle 
states into which the given three-particle state n may 
decay. 

The connected average of the five-particle term in 
Eq. (HI) with the same contribution from G'oo(e') leads 
to a pole at e' — £ = iT/5 in complete analogy with 
Eq. (pH). Here again r/5 is the smallest total width for 
the five particles with given total energy J2 — ^- The 
straightforward estimation of the five-particle contribu- 
tion to the return probability gives 

I^~(A3/r)'e-"/5 (16) 

Neither the overall numerical constant nor the "weak" 
~ 1/rt prefactor are shown here. Non-Cayley-tree con- 
tributions to P5 coming from the corrections to Gqq not 
included in the continued fraction Eq. (^) are F/A times 
smaller and decay also like e~^*^^ (or faster) at large time 
t. 

Generalisation of Eq. ( p^ ) for arbitrary number of final 
particles gives 

iW:'-(A3/r)"e-i^*/(2"+i) . (17) 

Adding one more particle-hole pair in the final state costs 
a small factor A3/r, but the time evolution of this contri- 
bution is governed by the smallest joint width for 2n + l 
particles with fixed total energy, which is F/ {2n+l). The 
total return probability is given by the sum over all many 
particle contributions ( p^ ) (starting from n — 0). How- 
ever, at any given moment of time one of the P2n+i dom- 
inates and all others may be neglected. Figure 1 shows 
the function In P(t) and illustrates how it is formed by 
different contributions from Eq. (p^). The two consec- 
utive values P2n+i coincides and the strength of decay 



is changed at Ft = {2n - l){2n + l)L/2. For large Vt 
one may also replace the piecewise linear InP by single 
smooth function (see Figure) such that 

(18) 

This formula should be compared with the usual decay 
exp(-F<) = cxp{-7r£2t/2g2A}. Eq. (^ is the main re- 
sult of this paper. Since the true many particle density of 
states is ^ exp(— 27r-\/e/6A) this decay may formally 
proceed until time t ~ {g^ /e)\/\n{e/ ^fglS). 

Within the interval 3L/2 < Tt < 15L/2 the second 
term in Eq. ( p^ dominates. Consider the distribution 
of P{t) in this region. The direct generalisation of the 
calculation described by Eqs. (|5|,p|,pl]) gives 

P{tY = n!P(7)" . (19) 

Here nl simply accounts for the number of ways how 
n Greens functions Goo(ei) may be contracted with n 
Gnn (e' ) . The distribution function corresponding to the 
eq. (y^ is evidently 

VF(P) = 1/Pexp{-P/P} . (20) 

Moreover, at any time only one contribution ([l^), cor- 
responding to fixed number of final particles, dominates 
in P{t). Therefore, Eq. ([l^) and consequently Eq. ( [20| ) 
arc also valid for any large time Ft > 3L/2. The strong 
fluctuations of P{t) in the asymptotics ( pO| ) are in con- 
trast with the usual decay of the return probability at 
rt < 3i/2, where one has varP(t) < P(t)^. 

It is seen from Eq. (|20| ) that the stretched asymptotic 
decay of a quasiparticle in a finite system described by 
Eq. (HI) does not require any special realization of the 
quantum dot. The interesting asymptotic effects due to 
rare fiuctuations of disorder were considered more than 
a decade ago in Ref . |l^ . (See also the recent discussion 
of the decay in case of strong interaction [p^ .) 

All the above results only correspond to the case of 
relatively high excitation energy. Below e ~ y^A the 
width expected from the golden rule become smaller 
then the level spacing for available three-particle states 
F < A3 (|). StiU at ^/g/lng < e/A < ^ quasiparti- 
cles are unstable . Not much is known about the time 
dependence of the return probability P(t) in this region. 
The decay now is not exponential at any time. One may 
convert the analysis of Ref. in order to get the time 
dependent series 

p2 /pi /AW" 

where the c„ are some (unknown) numerical coefficients. 
Unfortunately, Eq. ( pi] ) is valid only for 1 — P <C 1. This 
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result may be used for an estimation of the time after 
which the decay of the return probabihty starts, but it 
does not teach us about the functional form of P{t). 

To conclude, in this paper we have considered the de- 
cay of a single-particle excitation in a finite fermionic 
system with random two-particle interaction. Our results 
may be interpreted as follows. At short time t the decay 
into two-particles and one hole proceeds with the Breit- 
Wigncr width F, in accordance with the golden rule. At 
a certain time the probabilities to find the system in the 
initial or any of the discrete final three-particle states 
became of the same order of magnitude. After that the 
further decay (exponential on average) is governed by 
smaller width r/3, which is the minimal width of the cor- 
responding three-particle state. At these times the prob- 
ability P{t) to find the initial quasiparticle unperturbed 
fluctuates strongly, but we were able to find the distri- 
bution function W{P) describing these fluctuations (Eq. 
(pO|)). When the three-particle states become completely 
mixed with the five-particle states, the five-particle width 
come into play leading to P{t) ~ exp(— rt/5), and so on. 
This scenario with many stage decay takes place for the 
random interaction and an essentially irregular single- 
particle spectrum, which suppresses the interference be- 
tween different decaying channels. At sufficiently large 
time this piece-wise linear exponential decay may be de- 
scribed by the smoothed formula ( |l^ ) with ln(F) ^ — 
All these results are valid for high enough excitation en- 
ergy (e Vs^)- The behaviour of P(i) below the golden 
rule threshold, at e < y^A, requires a further investiga- 
tion. 
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FIG. 1. The logarithm of the return probability P as a 
function of Ft, both measured in units of large logarithm 
L — ln(r/A3) (thick polygonal line). The thin dashed lines 
are the individual contributions due to the three-particle de- 
cay {—Tt) and due to the correlations in the 2n -|- 1-particle 
final state {—Tt/(2n + 1)). The thick dashed line is the 
smoothed asymptotics P ~ exp{—\/2Lrt}. 
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